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1. Introduction 

Metrology is a prime issue both from theoretical and practical reasons. Precise 
measurements are crucial in physics since they constitute the link between the theory 
and nature, so that accurate measurements can promote or reject a theory. 

The key contribution of the quantum theory to metrology is that quantum 
fluctuations limit the resolution. To be more precise, quantum limits emerge when 
we impose constraints. The most popular restriction is to consider a fixed finite mean 
number of particles, although time limitations might also be considered pQ. 

The existence of quantum metrology limits imposed by finite resources is very clear 
if we consider systems with bounded number of particles [IH3] . In linear schemes this is 
the Heisenberg bound that states that the uncertainty is lower bounded by the inverse of 
the total number of particles employed in the measurement. However, most states have 
unlimited number of particles even if their mean number is finite. In such a case, despite 
that it may be expected a similar scaling law with the total number of particles as in 
the finite case, to give a ultimate bound to the sensitivity turns out to be considerably 
more intricate [4HT6]. 

In this work we provide evidence that the Heisenberg bound can be beaten with a 
fixed and finite mean number of particles. Sub-Heisenberg resolution arises because 
of the combination of two effects. Firstly, performance estimators (such as Fisher 
information [TTt[T8]) may be nonlinear functions of the photon number, even if we are 
dealing exclusively with linear processes. Secondly, the use of probes in states including 
coherent superpositions with the vacuum shifts the number distribution to larger photon 
numbers but keeps fixed the mean number. 

2. Quantum limits in the detection of non-random signals 

Throughout the article we consider that the signal is encoded as a shift in the phase of a 
quantum harmonic oscillator; a suitable practical implementation may be a single-mode 
electromagnetic field. Furthermore we consider non-random signals, i.e. the phase-shift 
is assumed to be a unknown but deterministic signal [T9H21] which is not subject to 
statistical fluctuations. Rather, it is the estimation process itself what is the source 
of randomness because the statistical nature of the quantum mechanics. For instance, 
parameters such as the mass of a particle or the amplitude of a gravitational wave (on 
a classical theory of gravitation) are commonly assumed to be deterministic signals. 
As a measurement of uncertainty we shall consider the mean squared error, 

(A0) 2 = J^dxp(x\<P)i4>(x)-(4>)] 2 , (1) 

here x denotes the possible outcomes of the performed measurement which takes value 
in the measurable set I. The outcomes x follow the distribution p(x\<f>) that depends on 
the true value of the phase-shift <fi. The function <f>{x) is an estimator of the phase-shift 
with expectation value (0) = J 7 dxp(x\(f))<f)(x) . For an unbiased estimator (0) = <fi. 
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Note that the sensitivity of a measurement may depend on the unknown value of 
the signal A0 = (A0)^. In fact, this is quite natural in real experiments (see [22] for 
example). Since quantum fluctuations are typically relevant only for small signals, we 
shall focus here on the behaviour of (A0)^ for small phase shifts <C 1. Moreover we 
are primarily concerned with the resolution for a specific value of the nonrandom signal, 
say = 0. This might be the case of the measurement of the photon mass M p h so that 
oc M ph [23]. Accordingly we care mainly for the uncertainty (A(f>)^ =0 at that specific 
point = disregarding its behavior at other signal values. 

In this regard, it is important to note the difference between (A0)^ and the averaged 
mean squared error, which is typical in Bayesian estimation approaches. In those 
situations, the signal is also considered to be a random parameter, and the uncertainty 
may be measured by the formula 



dxp(<p)p(x\<f))[<f)(x) — 0] = / dcf) I dxp(x, 4>)[4>(x) — 0] , (2) 

/ J-n Jl 

where p(4>) is the so-called prior probability. It may represent somehow the state of 
ignorance about prior to the experiment [21]. Note that for an unbiased estimator 
Eq. becomes 

50 2 = £ d<f>jp(x,<f>Mx) - (0)] 2 = ((A0) 2 >, (3) 

which coincides with the mean squared error flTJ) averaged with the prior probability. 

In recent papers [121416] several bounds have been established for the averaged error 
50 which behave like k/N^. Here k is a constant of the order of one and is the mean 
value of the total number of photons employed in the estimating procedure. 

However, as already mentioned, in this work we are interested in the non-averaged 
uncertainty A0. Previously [7] suggested that A0 is always larger or of the same order 
as l/iVr- In [9] it was proved that, within a two-field-mode context, A0 is lower 
bounded exactly by 1/iVx for states and measurements without coherences between 
different subspaces of fixed total photon number. That, in one-mode approximation, 
corresponds to l/(2iV T ). We shall refer to 1/ (2iV T ) as the "Heisenberg bound". For the 
general case [9] only argues that the Heisenberg bound strictly applies in the asymptotic 
limit when the number of repetitions becomes infinity. So A0 could still be lower than 
l/(2iVx) for finite number of measurements, as explicitly recognized in [9]. This is 
examined in Sections 3 and 4. 

For unbiased estimators, a very well-known inequality is the Cramer-Rao bound 
p/FlllE], which is formulated for (A0) 2 , 

^ y - - ^mr (4) 

where F(<ft) is the Fisher information of the measurement, m the number of repetitions 
of it, and -Fq(0) the quantum Fisher information which corresponds to the best possible 
choice of the measurement [18]. If such a measurement is good enough, this bound 
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can be achieved in the asymptotic limit m — > oo by using the maximum likelihood 
estimator [IHHZE]. 

For linear schemes, Fq(<P) is of the form an 2 + j3n, where usually a, (3 do not 
depend on the probe state. Hence (jlj) asserts that the optimum mean squared error must 
decrease as 1/n, which is in accordance to the result of [9]. Thus, we refer to A(f> oc 1/n 
(i.e. A<p oc 1/iVx) as "Heisenberg scaling". Sensitivities below the Heisenberg bound 
1/(2./Vt) may be compatible with the Heisenberg scaling, but the implication in the 
reverse sense is not true, so that if the Heisenberg scaling is beaten then the Heisenberg 
bound is certainly surpassed at least for some N T large enough. 

In principle, to break the Heisenberg scaling we may resort to non-linear schemes 
[3||24"] where Fq((f)) is polynomial of higher order in n. However, it might be broken also 
in linear schemes if a and/or /3 in Fq ~ an 2 + (3n depend on the mean number n of the 
probe state. We shall examine this possibility in Section 4. 

3. Probe states and signal transformation 

The structure of any signal-detection process, classical or quantum, is quite universal. 
A probe experiences a signal- dependent transformation. The change of the probe state 
is monitored by a measurement whose outputs serve to infer the value of the signal. 

Next we examine the construction of a probe state with the aim to beat the 
Heisenberg bound 1/(2 At) • 

3.1. Squeezed states 

As a first example let us consider the familiar choice |25j of a quadrature squeezed 
state [26], 

\O = D(y)S(r)\0), (5) 

with 

D(y) =exp(iyX/2) } S(r) = exp [ir (XY + YX) /4] , (6) 
where X, Y are the quadratures of a single-mode field of complex amplitude a, 

X = a) + a, Y = i (a f - a) , (7) 
with [X, Y] = 2i. Moreover we have 

<X) C = 0, {Y)t = y, (8) 

with 

v 2 

n ? = ( a t a ) s = t- + s inh 2 r, (AX)J = exp(-2r). (9) 

For definiteness, throughout we will assume y > 1 and r > 1. In addition we take 
a precisely equal splitting of the photons between the coherent and squeezed parts, 
y 2 /4 = sinh 2 r. Despite this is not the best distribution of the energy resources [27], it 
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does not qualitatively affect the results and simplifies subsequent computations. The 
following relations hold approximately 

1 



-2 2r 



y 



-2 ~ 2n € > 1. 



(10) 



(AX)J 

As usual, in linear schemes the phase shift to be detected cj> is generated by the 
number operator a^a so that the transformed probe state is exp(— i</>a'a)|£). In this 
situation the quantum Fisher information is four times the variance of a^a in the state 
If) B2], and yields 

F Qie = 4(Aa t a) 2 ~6n 2 . (11) 

In order to discuss whether the above result for the quantum Fisher information 
allows for an improvement on the Heisenberg bound, it may be convenient to consider 
a practical scheme, even if this is not fully optimal. To this end let us analyze 
the homodyne detection measuring the quadrature X in the transformed probe state 
exp(— i(jxrfa)\£). In this case the Cramer-Rao lower bound reads [17] 



~\ 2 1 
A0) > 



(12) 



J <t> mF((j)) ' 

where m is the number of repetitions of the measurement and the Fisher information 
reads 

>n 2 



dx- 



dp(x\(j)) 



p(x\4>) 

For a small enough signal cj> <C 1 we may approximate 

1 



(x\e 



10 



(2tt)V4 > /(AX) 



exp < i 



where 



and 



g (x) = - [x 2 



(AX)\ 



+ y- 



[x — x((p)]' 
4(AX)2 



x 



(AX); 



(13) 



,(14) 



(15) 



(16) 



x{4>) = ysincj) ~ ycj). 

Note that the change in (Y (((>)) = y + 0(4> 2 ) and AX(<p) = AX + 0(4> 2 ) is of second 
order in 0, because (X) = and (XY + YX) = 0. 

Under this assumption the statistics of the measurement is a Gaussian function 
centered in x(<p) = ycj), 



Pd x \ 



\(x\e-^ ta \0? 



cxp 



[x — x(cj>)f 
2(AX)2 



(17) 



leading to 



~\ 2 
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Since the estimation of <p is equivalent to the estimation of the mean of a Gaussian 
distribution, the maximum likelihood estimator reaches the above Cramer-Rao bound 
for any m [19-21]. Thus, the best that can be done with a squeezed probe is the 
single-trial sensitivity, m = 1, Nt = mn = n so that ( A0) ~ l/(2iVr), which is the 
Heisenberg bound. Therefore, to obtain sub-Heisenberg resolution, we must resort to 
another state. 



3.2. Superpositions of vacuum and Squeezed states 

Instead of a squeezed state, consider the probe in the superposition 

\j})=v\0) + v\O, (19) 

of the vacuum |0) and a squeezed coherent state |£) in the same terms as in the previous 
section. The parameters fi, v are assumed to be real for simplicity. Under the conditions 
( FlOl) the states |0) and |£) are approximately orthogonal 

(0|O = V / 2(AX) € exp [-f(AX)j/4] oc -L « 1, (20) 

so that fi 2 + v 2 ~ 1. Furthermore, we shall consider v -C 1, jJ, — 1. 
The mean number of photons of the state ffl9l) reads 

n = (^|a t a|V>) = (21) 

Since our purpose is to investigate quantum limits for fixed finite resources we will 
suppose that n is fixed, so that the number of photons in the squeezed state |£) depends 
on v as 

n, = -J. (22) 



In this case, considering 3> 1 and 

-V Tif — V Tit — 

2 ? i 2v 



(Aa^a) 2 c 5 -u 2 nl-^n 2 c% (23) 



and 



10f) 2 ~ v 2 

F ° = % ^ >- To^- < 24 > 

These results suggest that is able to estimate phase changes with an uncertainty far 
below the Heisenberg bound, as we have taken v <C 1. However, in order to confirm 
this, we need to find (if they exist) a good enough measurement scheme and estimator 
function. 

Let us consider again the homodyne detection measuring the quadrature X in the 
transformed probe state exp(— iifia 1 a) \ ip). Up to first order in <fr and v the statistics 
yields 

x 2 (x — y<j))' 



x 2 \ 2fiu _ [yx 



p(x\(f)) ~ fx exp I — — H ===== cos — — 4>g{x 



r*K V ~* V 2 ; ' ^(AX^ — L 2 ™~'J 6XP 1 4 4(AX)f 
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so that 



where g(x) is given by ( TT5"j) . 

This approximation can be checked by expanding p(x\4>) in power series of and 
comparing with the result of exp(— i<fia) a)\ip)\ 2 ~ — i^a^al^) | 2 at first order in 
v. Here for the vacuum we have 

( X |e-^|0> = (x|0) = ^-L 7I exp(-x 2 /4), (26) 

while Eq. ( Ti4l holds for the squeezed component. 

For the computation of the Fisher information, in the denominator in Eq. ffl3|) 
we can safely approximate p(x\4>) ~ (x|0) 2 . Taking into account relations (11 01) . i.e., 
y ~ 1/(AX)^ ~ a/2~^ ^> 1, we get that, after a long but straightforward calculation, 
the leading term in is 

F ~ 4z/ 2 n 2 ~ 4—, (27) 

This scales as the quantum Fisher information f )24|) . so that the conclusion is the same 
as above. 

For the sake of clarity, in Fig. 1 we have plotted p{x\4>) for n = 25, = 0, and 
v = 5 x 10~ 2 . It can be appreciated that the squeezed state contributes with a very 
narrow central Gaussian peak in a more uniform background provided by the vacuum. 
For small the effect of the phase shift becomes a displacement of the small peak 
proportional to 0. This peak is the cause of the increased Fisher information depending 
on v. However note that the Fisher information is only a lower bound for the resolution, 
which is only attainable for Gaussian distributions. Since the complete distribution (J25J) 
is manifestly nonGaussian, the lower bound cannot be achieved. Hence, the performance 
of the bound can only be tested by actual calculation with a specific estimator. That 
will be the aim of the next section. 

On the other hand, it is interesting to estimate how likely is the outcome x 
to fall below the tiny peak. Since the probability at x = is p(x = 0) ~ 
(ji 2 + 2fiu/ (AX)^j /V^tt, while its width is proportional to (AX)^, we obtain the 

probability below the small peak to be p(x = 0)(AX)^ ~ u/y^E+ z/ 3 / 2 /-^ 1 / 4 ~ u/y/E, 
where this approximation holds for v <^ 1. 

4. Results and Discussion 

In order to check that the state (fT§|) can actually beat the Heisenberg bound we have 
to choose an appropriate estimator. To that aim we consider the maximum likelihood 
estimator, = 0ml, such that 

£(0ml|zi,x 2 , • • • , x m) = max £(0| £1,0:2, • • • ,x m ), (29) 
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Figure 1. Probability distribution of the quadrature X in the probe state (fT9| as a 
function of x for n = 25 and v = 0.05. The inset represents in more detail the central 
peak. 

where the likelihood function is 

in 

£((f)\xx,X2, . . . ,x m ) = [[p(xi\(f>). (30) 

i=i 

The maximum likelihood estimator is unbiased for any value of <j) and it is asymptotically 
efficient [T9H2T] . i.e. it reaches the Cramer- Rao bound for a large number of 
measurements m — > oo. 

However for large m, Monte Carlo simulations seem to indicate that the sensitivity 
is worse than the Heisenberg bound. That basically happens because the uncertainty 
decreases with m as l/y/m whereas the Heisenberg bound does it as 1/m. Thus we 
restrict our study to the best situation possible where m — 1. The table 1 shows some 
of the uncertainties obtained by Monte Carlo simulations for different mean number of 
photons. 

As it can be seen, the probe state f|T9l) provides sensitivities considerably smaller 
than the Heisenberg bound, but compatible with the Cramer- Rao bound f )28|) . For 
instance, for m = n = 1, v = 0.05, it reads (A0)^ = o > 0.025. On the other hand, 
in Fig. 2 we have plotted the evolution of the uncertainty (A0)0 =o with the number 
of points taken to simulate the experiment with n = 1 by the Monte Carlo procedure. 
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n 


d> 




1/(2JV T ) 


1 


1.60 x 10" 4 


0.035 


0.500 


2 


-9.50 x 10~ 5 


0.025 


0.250 


3 


8.79 x 10" 5 


0.020 


0.167 


4 


-5.43 x 10~ 5 


0.017 


0.125 


5 


4.91 x 10~ 5 


0.016 


0.100 



Table 1. Results of the Monte Carlo simulation for (ft — and different mean number 
of photons. It has been taken v = 0.05 and m = 1. The two last columns show clear 
violations of the Heisenberg bound l/(2iV T ). 
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Figure 2. Evolution of the uncertainty (A<ft)^, = o with the number of points of the 
Monte Carlo simulation for m = 1, v = 0.05 and h = 1. 



It clearly shows that statistical uncertainties are small enough to reject them as the 
possible reason for the violation of the Heisenberg bound. 

Similar uncertainties are obtained up to the fourth decimal place for any value of 
(f) inside of (A^)^, =0 . Of course, for the range of parameters we have considered, it is 
not possible to distinguish between such a finite (ft and = 0. However this state is just 
intended to be an example of "proof of principle" that the Heisenberg bound may be 
suppressed. 
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Figure 3. Logarithmic plot of the sensitivity (A(j))^ = Q with m = 1 and v — 0.05/n as 
a function of n, for some values between 1 and 5. The linear fitting (red line) leads to 
a sub-Heisenberg scaling (A0) 0=o = (0.0354 ± 0.0003)/7V^ 4976±0 0098 . 



Interestingly, we may also explore whether the probe state (JT9l) can beat the 
Heisenberg scaling. At difference with previous approaches, in our case the uncertainty 
depends on an extra free parameter v that depends on the probe state. Thus, we are 
free to consider that v might depend on n for some probe states. Note that provided 
the mean number of squeezed photons is given by Eq. (I22I) . = n/u^n) 2 , the total 
number of photons n is not affected by v(n) despite of its dependence with n. Fig. 3 
shows the results of a Monte Carlo simulation for v = 0.05/n. The linear fitting of the 
data provides a scaling law of the form (A0)^ =o ~ i/jv^- 4976 * - 0098 , which improves the 
Heisenberg scaling (A0)^ =o ~ l/N^. 

It is remarkable that in [4] another scheme which seemed to beat the Heisenberg 
bound was proposed. Incidentally our Fig. 1 presents some similarities with the figure 
in [1]. However, in contrast with our case, Shapiro et al. jl] focus on the results of a 
covariant measurement, where the uncertainty is the same for any value of <fi. For that 
case, Braunstein et al. [5] showed that the sensitivity neither improves the Heisenberg 
bound nor the Heisenberg scaling. This is actually consistent with the results in [T2HT6] , 
where the averaged mean squared error (T5]) is analysed. Our approach is different from jl] 
and [5] in the sense that we do not look for a 0— independent sensitivity, but a large 
sensitivity at a very small interval around one point (namely = 0). Because of this 
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we can resort to homodyne measurements which are not covariant. 

Thus, we want to emphasize that the improved resolution holds only for <p close 
enough to <fi = 0. Nevertheless, the limitation of improved resolution to small intervals of 
signal values is quite frequent in practice. For example in [22J improved sensitivity holds 
just for a signal interval 8<p of the same order of the uncertainty 5<p — A0. Therefore, 
it appears that the condition A<p <C 5<p is not mandatory in real experiments. 

5. Conclusion 

Quantum metrology is framed by some implicit understandings that may limit its 
development. In this work we have shown that sub-Heisenberg resolution is possible for 
local estimators. Concretely we give an example where homodyne measurement provides 
uncertainties smaller than 1/(2JVt) for signals close enough to <p = 0. This has to 
be contrasted with previous results involving covariant measurements schemes [HE], 
asymptotic number of repetitions |9j, or the recent rigorous proof of sensitivity bounds 
for averaged mean squared errors in Bayesian estimation [T2T4T6] . 

The reason for this improvement lies in the nonlinear behavior of performance 
estimators with the number of photons. This leads to resolution improvement at fixed 
mean number of photons because of the shift in the photon-number statistics caused by 
the coherent superposition with the vacuum. 

Local measurement and uncertainties may be especially useful if the experimenter 
knows a priori that the signal is very close to a certain value. However our aim here 
has been just to provide evidence that resolutions beyond the Heisenberg bound are 
possible. We hope that may disclose a new perspective on quantum metrology avoiding 
previously assumed performance limits. 
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